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ABSTRACT 

More than a dozen blazars are known to be emitters of multi-TeV gamma rays, often with strong and rapid 
flaring activity. By interacting with photons of the cosmic microwave and infrared backgrounds, these gamma 
rays inevitably produce electron-positron pairs, which in turn radiate secondary inverse Compton gamma rays 
in the GeV-TeV range with a characteristic time delay that depends on the properties of the intergalactic mag- 
netic field (IGMF). For sufficiently weak IGMF, such "pair echo" emission may be detectable by the Gamma- 
ray Large Area Space Telescope (GLAST), providing valuable information on the IGMF. We perform detailed 
calculations of the time-dependent spectra of pair echos from flaring TeV blazars such as Mrk 501 and PKS 
2155-304, taking proper account of the echo geometry and other crucial effects. In some cases, the presence 
of a weak but non-zero IGMF may enhance the detectability of echos. We discuss the quantitative constraints 
that can be imposed on the IGMF from GLAST observations, including the case of non-detections. 

Subject headings: BL Lacertae objects: general — gamma rays: theory — magnetic fields — radiation mech- 
anisms: non-thermal 



1. INTRODUCTION 

Magnetic fields are observed to be common in the struc- 
tured regions of the universe, such as galaxies and clusters 
of galaxies. To interpret their microgauss field strengths, 
the pre-collapse seed fields in the intergalactic medium 
required by dynam o theories may be of order 10~ 20 G 
dKulsrud et al. 19 97). Various theoretical possibilities have 
been suggested for the origin of such large-scale, inter- 
galactic magnetic fields (IGMFs). For example, primor- 
dial magnetic fields of ~ 10~ 20 G may b e generated dur - 
ing the cosmological QCD phase transition (Sigl et al. 1997b . 
They can also be produced during the inflationary epoch if 
conf ormal invariance of elect romagnetic interactions is bro- 
ken dTurner & Widrow 19 88). Density fluctuations before 
cosmic recombination inevitably give r ise to weak IGMFs 
(Taka hashi et al. 2005Ulchiki et al. 20061) . The Biermann bat- 
tery mechanism ((Gnedin et al. 20001) or radiation drag effects 
(lLanger et al. 20051) at cosmic reionization fronts can induce 
fields of 10" 20 - 10" 16 G. At low redshift, such weak IGMFs 
are expected to survive in intergalactic void regions, as they 
can remain uncontaminated by ast rophysical sources such a s 
galactic winds or quasar outflows (Furl anetto & Loeb 2001). 
Measurements of IGMFs would be crucial for understanding 
the origin of galactic magnetic fields. 

One of the best-known tools to probe magnetic fields is 
Faraday rotation measurements, from which an upper limit 
of ~ 10~ 9 G has been der ived for the IGM F assuming cor- 
relation lengths ~ 1 Mpc (Kronberg 19941). A different ap- 
proach suitable for probing very weak IGMFs is to uti- 
lize "pair echo" emission from transient very-high-energy 
(VHE) gamma- ray emitters such as blazars and gamma-ray 
bursts (GRBs) dPlaga 1995b . Multi-TeV photons from dis- 
tant sources are attenuated by 77 pair-production interactions 
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with the cosmic microwave background (CMB) and cosmic 
infrared background (CIB). The electron-positron pairs then 
up-scatter CMB and CIB photons by the inverse-Compton 
(IC) process to produce secondary gamma rays, whose flux 
depends strongly on the IGMF. If it is stronger than ~ 10~ 12 G, 
the formation of a very extended and ne arly isotropic pair halo 
is unavoidable ( Aharonia n et al. 19941) . On the other hand, if 
the IGMF is sufficiently weak, most of the secondary gamma 
rays will come from the direction of the source with some 
temporal and spatial spreading as pair echos, containing valu- 
able information on the IGMF. Blazars are promising sources 
for this purpose, since they are observed to exhibit strong 
flares at multi-TeV energies with flux variations by a factor 
of > 10 over time scales of < 1 hr to months. Pair echo emis- 
sion is typically expected at GeV energies, appropriate for the 
recently launched Gamma-ray Large Area Space Telescope 
(GLAST) satellite. 

Here we reconsider pair echo emission from blazars by 
exploiting the formulation recently developed by Ichiki et 
al. (2008; see also Takahashi et al. 2008), which al- 
lows more satisfactory calculations of the time-dependent 
echo spectra c ompared to previous works dDai et al. 20021 : 
iFan et al.T0 04). Since the pair echo emission can persist af- 
ter the primary flare decays, it may be observable unless it is 
hidden by some quiescent emission. Furthermore, constraints 
on the IGMF are possible even in the case of GLAST non- 
detections. We demonstrate this using as examples the past 
flares from Mrk 501 in 2005 and from PKS 2155-304 in 2006. 
Further details of our methods and results will be described in 
a subsequent paper (Takahashi et al, in preparation). 

2. EMISSION PROPERTIES 

So far, 23 AGNs have been detected at energies > 0.1 TeV 
4 (Wagner 2008). Most of them belong to the high-frequency- 
peaked BL Lac subclass of blazars, characterized by spectral 
energy distributions with two maxima occurring in the X-ray 
and TeV gamma ray bands. The standard blazar model com- 
prises supermassive black holes ejecting relativistic jets close 

4 http://www.mppmu.mpg.de/ rwagner/sources/ 
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to the line of sight. Owing to relativistic beaming, blazars 
exhibit rapid flux variations on timescales down to a few min- 
utes, with strong TeV flares being the most extreme events. 

The 77 optical depth of the CIB depends on gamma- 
ray energy, source redshift and CIB intensity. For the 
CIB, we adopt here the "best-fit" and "low-IR" models of 
Kneiske et al. (20 02, 2004) (see also iPrimack et al. 20051: 
IStecker et al. 20061) . Note that recent observations of TeV 
blazars may point to a CIB resembling the low-IR model, 
close to the lower limits from galaxy count data (Aharonian et 
al. 2006; Albert et al. 2008; see however, Stecker et al. 2007). 

2.1. Pair Echo Emission 

The total fluence of pair echo emission is determined by the 
77 optical depth of the CIB and does not depend on the IGMF. 
Primary photons with energy E 1 are converted to electron- 
positron pairs with Lorentz factor 7,, w 10 (£' 7 /lTeV)(l+z) 
in the local cosmological rest frame, which then up-scatter 
CMB and CIB photons. CMB photons are boosted to ener- 
gies ~ 2.82fe^ MB7e 2 /(l +z) » 0.63(£ 7 /l TeV) 2 (l +z) 2 GeV, 
where 7^ MB ~ 2.73(1 +z) K is the local CMB temperature. 
To evaluate the pair echo flux, we must consider various time 
scales involved in the process, such as the flare duration, angu- 
lar spr eading time, and the delay time d u e to magnetic deflec - 
tions dDai et al. 20021: iFan et al. 20041 iMurase et al. 20071) . 
These can be estimated as follows. 

The angular spreading time is Af ang w (1 + 

z)(A; c + A 77 )/2 7e 2 c, where A 77 « (0 .26a T n' cm y l « 

20Mpc(«Q B /0.1 cm" 3 ) 1 is the local 77 mean free path 
in terms of the local CIB photon density n' cm , and 

AJ C = 3m e c 2 /(4a T U CMB j e ) w 690kpc(7 e /10 6 ) _1 (l + z)" 4 
is the local IC cooling length in term of the local CMB 
energy density U CMB . At the energies of our inter- 
est, A 77 > A{ c so that Af ang rj (1 + z)A 77 /2 7c 2 c « 

960s( 7f ,/10 6 )" 2 (n^ IB /0.1cm- 3 )" 1 (l + z). For suf- 

ficiently small deflections in weak IGMFs with 
present-day amplitude Big = B' ic (l+z)~ 2 and coher- 
ence length A co h = A£ oh (l + z), the magnetic deflec- 

1 li 

tion angle is 9 B = max[A{ c /r L ,(A; c A^ oh ) /r L ], where 
fL = "feiiieC 2 /eB[ G is the Larmor radius of the electrons or 
positrons. The delay time due to magnetic deflections is 
At B ~ (1 +z)(A( c + X' )(0l/2c). Note that prior to Ichiki 
et al. (2008), the A 77 term here had been neglected and 
At B underestimated by up to 2-3 orders of magnitude. For 
coherent magnetic fields with X' coh > AJ C , we have At B w 

max[6.1 x 10 3 s(% / 10 6 )' 5 (B 1G /lO' 20 Gf (I +zT\ 1. 6 x 

10 5 s( 7 ,/10 6 )" 4 («ciB/0.1cm- 3 )" 1 (BiG/10- 20 G) 2 (l+zr 3 ]. 
Implicit in the above discussion is that both 1/7^ and 9 B do 
not exceed 9j, the opening angle of the AGN jet; otherwise a 
significant fraction of photons or pairs will be deflected out 
of the line of sight and the echo greatly diminished. 

Together with the flare duration T, the pair echo delay 
time can be estimated by Af = max[Af ang , At B ,T]. If 
At B dominates, pair echos can serve as effective probes 
of weak IGMFs. For flaring blazars with T ~ day, 
Af ang is irrelevant, and At B > T at echo emission energies < 

0.86GeV(B IG /10- 20 G)(r/l day)~ 1/2 (n CIB /0.1 cnr 3 )~ 1/2 (l +z)~ 
for the case of coherent IGMFs. Such estimates were 
the basis of p r evious evalua tions of the pair echo flux 
dDai et al. 20021: IFan et al. 20041) . but explicit descriptions of 



the time-dependent sp ectra were not po ssible without some 
ad hoc modifications (Murase 'et al. 2007h . In constrast, the 
formulation of Ichiki et al. (2008) enables us to calculate the 
time-dependent spectra in a more satisfactory manner, partic- 
ularly at late times, accounting properly for the geometry of 
the pair echo process. 

2.2. Primary Emission Spectrum 

First we consider the archetypal flaring blazar Mrk 501. 
Strong flares at energies up to 20 TeV were observed in 1997 
by HEGRA (e.g., Katarzynski et al. 2002). Similar strong 
flares wer e recently observed by MAGIC from May through 
July 2005 (lAlbert et al. 20061) . during which the flux varied by 
an order of magnitude, and intranight variability was observed 
with flux-doubling times down to 2 minutes on the nights of 
June 30 and July 9. We focus on the strong flare of June 30. 

The second example is PKS 2155 -304, the brightest VHE 
blazar in the southern hemisphere dAharonian et al. 2 005a; 
I Aharonian et al. 2005bT) . During July 2006, the average VHE 
flux observed b y HESS was more than 10 times its quiescent 
value in 2003 (Aharonia n et al. 20071) . In particular, an ex- 
tremely strong flare was observed on July 28, ~ 50 times 
brighter than the quiescent level, which we use as a tem- 
plate for calculating the pair echo emission. Only small spec- 
tral differences were found between the flaring and quiescent 
states, as opposed to other blazars t hat often reveal large sp ec- 
tral changes at different flux levels dAharonian et al. 2 007). 

From a theoretical viewpoint, the intrinsic maximum en- 
ergy £™ ax should reflect either the maximum energy of ac- 
celerated electrons or protons, or a cutoff due to internal 77 
absorption. Although the true value of E™ ax is not yet known, 
the spectrum of Mrk 501 was observed to extend at least to 
~ 20 TeV, so here we take £ , ™ ax = 20 TeV as a reasonable as- 
sumption. Concerning the emission mechanism, both leptonic 
and hadronic models have been proposed (see, e.g. Sikora & 
Madejski 2001; Btittcher 2006 for reviews). The synchrotron 
self Compton (SSC) model is one of the most popular lep- 
tonic models (e.g., Maraschi et al. 1992; Inoue & Takahara 
1996). Another frequently discussed leptonic model is the 
external Compton model, where electrons up-scatter exter- 
nal photons originating outside the jet. In BL Lac objects, 
the tightly correlated variability in the X-ray and TeV bands 
(e.g., Katarzynski et al. 2005) and the lack of strong emission 
lines indicate a minor role for external photons and favor the 
SSC model. In hadronic models involving protons accelerated 
to ultra-high energies, the high-energy spectra are attributed 
to synchrotron radiation from either the protons themselves 
(Aharonian 2000), or secondary electron-positron pairs gen- 
erated in photohadronic interactions (Mannheim 1993). The 
hadronic models are challenged by the observed X-ray -TeV 
correlation and rapid gamma-ray variability, but they have not 
been entirely ruled out. 

3. RESULTS 

Fig. 1 shows the time-dependent pair echo spectra for the 
2005 flare of Mrk 501 at different times t after the onset of 
the flare, whose flux is assumed to decay as exp(-f/r) on 
a timescale T' = T/(l+z) = 0.5 day. We take A coh = 0.1 
kpc and different values of Big- The sub-TeV primary spec- 
trum of Mrk 501 is relatively hard with pho ton index a < 
0.4 dAharonian et al. 20021: lAlbert et al. 20061) . Extrapolating 
this to GeV, the pair echo should be visible relative to the pri- 
mary flare. It should not be masked by the quiescent GeV 
emission either, if the latter is estimated by a one-zone SSC 
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FIG. 1. — Pair echo spectra for the 2005 flare of Mrk 501, plotted at 
t = T = 0.5 day (thin) and t = 3T = 1.5 day (thick), for the cases of B IG = lO" 20 
G (solid) and 10~ 18 G (dot-dashed) with A co h = 0.1 kpc. Also shown are the 
observed primary spectrum (thick dashed) and intrinsic primary spectrum for 
the low-IR CIB model (thin dashed) at t = 0, described by linear extrapola- 
tion at < 200 GeV. The quiescent emission is represented by an SSC model 
(dotted). The GLAST sensitivity for integration time t = 3T = 1.5 day is 
overlayed. 
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FIG. 2.— Pair echo spectra for the 2006 flare of PKS 2155-304, plotted at 
t = T = 0.5 day (thin) and t = 3T = 1.5 day (thick), for the cases of A coh = 0. 1 
kpc (solid) and A co h = 1 Mpc (dot-dashed), with Biq = 10~ 20 G. Also shown 
are the observed primary spectrum (thick dashed) and the intrinsic primary 
spectrum for the low-IR CIB model (thin dashed) at t = 0, described by a 
hadronic model at < 200 GeV. The quiescent emission is represented by an 
SSC (thin dotted) or a hadronic (thick dotted) model. The GLAST sensitivity 
for integration time t = 3T = 1.5 day is overlayed. 



TABLE 1 

Lower bounds on the IGMF in the case of GLAST 

NON-DETECTIONS OF PAIR ECHOS 



model consistent with the "low flux" TeV data of 2005 as in 
Fig. 21 of Albert et al. (2007). Compared with the GLAST 
sensitivity 5 , the detection prospects for pair echos seem rea- 
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FIG. 3. — Primary and pair echo light curves for a hypothetical strong flare 
of PKS 2155-304 compared with the GLAST sensitivity at 1 GeV (thick) 
and 10 GeV (thin), for the case of Biq = 10~ 20 G (solid) and B IG = G (dot- 
dashed) with A co h = 1 Mpc. The best fit CIB model is employed. Also shown 
is the total primary light curve for both flare and quiescent emission assuming 
an SSC model. 

sonable as long as the IGMF is sufficiently weak (see also 
Dai et al. 2002). Detailed observations of its spectra and light 
curve will provide valuable information on the IGMF (Taka- 
hashi et al., in preparation). 

Note, however, that Mrk 501 has been previously detected 
at GeV by EG RET during a mu ltiwavelength campaign in 
March 1996 dKataoka et al. 19991) . The TeV spectrum was 
much harder than in 1997, and the entire GeV- TeV spectrum 
was incompatible with the simplest, one-zone SSC model. If 
this GeV emission corresponds to a persistent, quiescent com- 
ponent of Mrk 501, it may obscure the pair echo emission. 
This question should be resolved soon by GLAST. 

Results for the 2006 flare of PKS 2155-304 are shown in 
Fig. 2, with assumptions for the flare similar to Mrk 501. 
Here we fix Big = 10~ 20 G and show the dependence on A CO h- 
The case of A co h = 0.1 kpc corresponds to tangled IGMFs 
(Acoh ^ Aic), while A co h = 1 Mpc implies coherent IGMFs 
(Acoh ^ Aic)- The quiescent GeV emission is estimated with 
both SSC and hadronic models for the average TeV flux in 
2003 as in Fig. 10 of Aharonian et al. (2005b). The flare 
spectrum of PKS 2155-304 above TeV is much softer tha n 
in Mrk 501 (lAharonian et al. 2005bl: lAharonian et al. 2007h . 
This may also be true b elow TeV if the SSC model is adopted 
(lAharonian et al. 20061) . in which case the pair echo may be 
hidden beneath the primary and quiescent emission. 

Generally speaking, stronger IGMFs dilute the echo emis- 
sion and make its observation more difficult. However, in 
some cases, the presence of the IGMF can enhance the de- 
tectability of the echo emission, at least in principle. As an 
illustrative example, we show the results for a hypothetical 
TeV flare of PKS 2155-304 that is a further 10 times stronger 
than in Fig. 2, using the best fit CIB model. Here the pair echo 
emission below ~ GeV outlasts the primary flare by virtue of 
the IGMF, enabling the echo to emerge after the flare subsides. 

Next we discuss the lower limits that can be imposed on 
the IGMF even when GLAST does not detect pair echos 
from TeV blazars. Such limits would be valid if the pri- 
mary flare and quiescent emission at GeV energies are low 
enough so that the case of = implies an excess of the 
echo flux dF l sec /dEj over the primary flux dF y ?r i/dE 7 . The 
non-detection of the echo emission can then be attributed to 
the effects of a finite IGMF, expressed as {dF 1 , sec / 'dE^) < 
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m?Lx[(dF 1 ^i/dE 1 ) 1 (dF 1 k m /dE 1 )\, where dF 1 ,\\ m /dE 1 is the 
GLAST sensitivity. Summarized in Table 1 are the constraints 
thus derived for the IGMF, which depend on the CIB as well 
as on assumptions for the primary emission. For PKS 2155- 
304, the limits are given only for the hadronic model, since the 
Big = echo flux is not expected to exceed the primary flux for 
the SSC model. More conservative but less model-dependent 
constraints may be deduced forblazars such as Mrk 501 with 
hard primary Te V spectra, where the primary GeV emission is 
expected to be less obstructive (Fig. 1). Note that these limits 
are for tangled IGMF that lead to lower echo fluxes at < 10 
GeV and hence more conservative limits compared to coher- 
ent fields. The same is true for low-IR models compared to 
best-fit models for the CIB. A more detailed account of the 
primary light curve should allow more realistic constraints. 

4. SUMMARY AND DISCUSSION 

We have evaluated the time-dependent spectra of secondary 
pair echo emission from TeV blazars and discussed the in- 
formation that can be derived for the IGMF, applying a re- 
cently developed formalism of pair echos that properly de- 
scribes their time evolution. The observational prospects 
are quite interesting for the recently launched GLAST mis- 
sion, and successful detections would open a new window 
on studies of cosmic magnetic fields. Even in the case 

of non-detections, lower limits on the IGMF of Z?igA 



1/2 



> 



-2K 



(10~ 19 - 10" 21 )GMpc 1/2 may be obtained, making use of suit- 
ably strong TeV flares with hard spectra. The existence of a 
weak but non-zero IGMF may also sometimes enhance rather 
than diminish the detectability of echo emission at late times. 

Further detailed calculations utilizing Monte Carlo meth- 
ods may be desirable for more robust predictions. For ex- 
ample, the effect of cooling of pairs during propagation in 



the IGMF can be moderately important, leading to fluxes 
small er by a factor of se veral compared to the results given 
here dMurase et al. 2007b . Also of concern are uncertainties 
in the CIB models, which can affect not only the pair echo flu- 
ence but also the timescales for angular spreading and mag- 
netic deflection delay at all redshifts. We must also beware 
of uncertainties in the intrinsic primary spectra including the 
value of £™ x . 

In addition to pair echos from flares, the quiescent emis- 
sion of TeV blazars may also contain useful information on 
the IGMF. For example, depending on the CIB, the spec- 
tra corrected for 77 absorption for some blazars including 
Mrk 501 point to a sharp pile-up at high energies, contra- 
dicting the expectations from conventional emission models. 
However, as long as Big ^ 10~ 18 G and the primary spec- 
tra has photon index a ~ 2 and £™ ax > 100 TeV, an inter- 
galactic cascade component may contribute to the quiescent 
TeV emission and co mpensate the effects of CIB absorption 
dAharonian et al. 20021) . The quies cent flux at < GeV could 
also be affected by such cascades dDai et al. 2002h . These is- 
sues will be elaborated on in a future publication (Takahashi 
et al., in preparation). 
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